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ABSTRACT 

GRB 080913 and GRB 090423 are the most distant gamma-ray bm'sts (GRBs) known to-date, with 
spectroscopicahy determined redshifts oi z — 6.7 and z = 8.1, respectively. The detection of bursts at 
this early epoch of the Universe significantly constrains the nature of GRBs and their progenitors. We 
perform population synthesis studies of the formation and evolution of early stars, and calculate the 
resulting formation rates of short and long-duration GRBs at high redshift. The peak of the GRB rate 
from Population II stars occurs at z ~ 7 for a model with efficient /fast mixing of metals, while it is 
found at z ~ 3 for an inefficient/slow metallicity evolution model. We show that in the redshift range 
6 ^ z ^ 10 essentially all GRBs originate from Population II stars, regardless of metallicity evolution 
model. These stars (having small, but non-zero metallicity) are the most likely progenitors for both long 
GRBs (coUapsars) and short GRBs (NS-NS or BH-NS mergers) at this epoch. Although the predicted 
intrinsic rates of long and short GRBs are similar at these high redshifts, observational selection effects 
lead to higher (factor of ^ 10) observed rates for long GRBs. We conclude that the two recently observed 
high-z GRB events are most likely long GRBs originating from Population II coUapsars. 

Subject headings: binaries: general — gamma rays: bursts — stars: formation 



1. INTRODUCTION 

The Swift satellite has recently discovered two high- 
redshift GRBs: GRB 080913 at redshift z = 6.7 (Schady et 
al. 2008), and GRB 090423 at redshift z = 8.1-8.3 (Tan- 
vir et al. 2009, Salvaterra et al. 2009). The gamma-ray 
properties of these bursts straddle the traditional dividing 
lines between short and long bursts, making a solid classi- 
fication of each burst difficult. For example, the observed 
burst duration of GRB 080913, 8 s, classifies it as a 
long GRB (Stamatikos et al. 2008). However, in the rest 
frame of the burst the duration is '--^ 1 s, suggesting a short 
GRB (Perez- Ramirez et al. 2008). It has been argued 
that, despite its relatively short duration (in the comov- 
ing frame), GRB 080913 nonetheless should be classified 
as a long GRB due to the specifics of the category defini- 
tion (which places the dividing line at 2 s in the observer 
frame, for a GRB sample with mostly unknown redshifts), 
and because other properties of this burst (e.g., the lag- 
luminosity and Amati (2006) relations) are consistent with 
those of long GRBs (Greiner et al. 2008). Similarly, GRB 
090423 had a burst duration of 9 s, and a rest-frame 
duration of ^ 1 s. Although many of the gamma-ray di- 
agnostics of this burst (Tgo, Epoak, hardness ratio) would 
classify it as "short" had the burst occurred at low red- 
shift, it has been argued that GRB 090423, as in the case 
of GRB 080913, is a "long" burst (Salvaterra et al. 2009; 
Zhang et al. 2009). 

We employ two simplifications of nomenclature. First, 
we retain the use of the "short" and "long" classifications 
for GRBs, despite growing evidence that the classification 
should be determined using a broader set of criteria beyond 



merely burst duration (Donaghy et al. 2006; Zhang et al. 
2007; Bloom et al. 2008). Second, we tie these two clas- 
sifications to specific progenitors using the simple theory 
correspondences outlined in Popham et al. (1999). Long 
GRBs are believed to be associated with the deaths of 
massive stars, and the associated formation of black holes 
(e.g., Woosley 1993; Galama et al. 1998; Hjorth et al. 
2003; Stanek et al. 2003; Woosley & Bloom 2006). Short 
GRBs, on the other hand, are thought to originate from 
the merger of compact objects, such as double neutron star 
(NS-NS) or black hole neutron star (BH-NS) binaries (e.g., 
Paczynski 1986; Eichler et al. 1989; Nakar 2007 and ref- 
erences therein). An alternative, unified model for GRBs 
has been proposed to simultaneously explain both long 
and short GRBs in terms of the same underlying physical 
mechanism, involving a black hole as central engine (e.g., 
RufRni et al. 2006; Dar & Rujula 2004). In addition, 
there exist models for both short and long GRBs that do 
not involve black holes (e.g., Usov 1992; King, Olsson, & 
Davies 2007; Cheng & Dai 1996). A preferable classifica- 
tion of bursts might focus on the progenitor classification 
of each bursts, which depends upon supernova association, 
host galaxy type, metallicity, and offset of the GRB with 
respect to the host galaxy (see Fryer et al. 2007 for a 
review). In this paper we study the progenitors directly, 
and identify coUapsars with long GRBs, and compact ob- 
ject mergers with short GRBs. 

Given their cosmological nature, and their bright af- 
terglows, GRBs are unique tools to probe cosmological 
parameters and the structure and composition of matter 
along the GRB sightlines (e.g., Lamb & Reichart 2000; 
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Lloyd-Ronning et al. 2002; Gou et al. 2004; Prochaska 
et al, 2007; Lamb 2007; Hartmann 2008; Savaglio et al. 
2008; Piro et al. 2008; Hartmann et al. 2009). This argu- 
ment gains force with the discovery of GRBs at very high 
redshifts; for example ORB 050904 at z — 6.29 (Cusumano 
et al. 2006; Kawai et al. 2006), and now GRB 080913 and 
GRB 090423 beyond the furthest quasars (e.g., CFHQS 
J2329-0301 at z = 6.4; Willot et al. 2007) and the fur- 
thest spectroscopically confirmed galaxies (galaxy lOK- 
1; J132359.8-f272456, at z = 6.96; lye et al. 2006). 
These high redshift objects provide a unique opportu- 
nity to study the Universe during the extended, inho- 
mogeneous epoch of reionization, tracing cosmic star for- 
mation back to the first generation of stars (Population 
III). GRB follow-up observations have the potential to 
reveal the fraction of neutral hydrogen and metal abun- 
dances as a function of redshift, possibly beyond z ~ 10, 
with future instrumentation on large aperture telescopes 
on the ground and in space. GRBs constrain star forma- 
tion models in the pristine Universe, and serve as unique 
light sources to pinpoint star formation activity even in 
the smallest proto-galaxies which would otherwise remain 
undetected. Rapid absorption spectroscopy of their af- 
terglows offers a truly exceptional opportunity for stud- 
ies of cosmic chemical evolution in an otherwise inaccessi- 
ble part of the cosmic baryon field. However, in order to 
utilize these distant GRBs as probes, prompt localization 
and redshift determination is crucial. Currently, the Swift 
satellite performs rapid response observations in the X-ray 
and optical-UV bands. Future missions, such as the Ener- 
getic X-ray Imaging Survey Telescope (EXIST; Grindlay 
2006), will provide combined wide-field X-ray and near- 
infrared imaging spectroscopy to study distant GRBs and 
associated black hole formation in the early Universe. 

In this study we attempt to characterize the progeni- 
tors of GRB 080913 and GRB 090423, and other potential 
GRBs at these high redshifts. As we shall see, this epoch 
is entirely dominated by Population II stars and we focus 
on GRBs produced only from these stars. We analyze the 
formation of stars at high redshift in §2. In §3 we use 
a self-consistent model to follow in detail the evolution 
of the most likely GRB progenitors: massive stars and 
NS-NS/BH-NS mergers. In §4 we adopt a cosmological 
model to estimate the GRB intrinsic (§ 5) and observed 
(§6) event rates at high redshift, explicitly identifying the 
likely progenitors of GRB 080913 and 090423. We con- 
clude in § 7. 

2. STELLAR POPULATIONS: MODEL/RESULTS 

Star Formation History. To derive rates for various stel- 
lar events (supernovae, mergers, etc.) one must start from 
a star formation rate history, and add assumptions about 
the cosmic evolution of basic properties such as the Initial 
Mass Function (IMF) and the fraction of stars in binaries. 
We use the analytic cosmic star formation rate (SFR) pre- 
scription provided by Strolger et al. (2004): 

sfr(t) = 10^a(i''e-*/'= -f de'^(*-*'')/'=) Mq yr^^ Gpc"^ (1) 

where t is the age of Universe (in Gyr) as measured in 
the global comoving (rest) frame, to is the present age of 
the Universe {to ~ 13.47 Gyr; see §4), and we adopt pa- 
rameters corresponding to the extinction-corrected model: 



a = 0.182, b = 1.26, c = 1.865, d = 0.071. The SFR 
rate density described above is in comoving units (space 
as well as time); for a non-evolving population, this rate 
is constant. 

Galaxy Mass Distribution. At redshift z < 4 we de- 
scribe the distribution of galaxy masses using a Schechter- 
type probability density function, calibrated to observa- 
tions from Fontana et al. (2006): 

$(Mgai, z) = $*(z) ln(10)ai+"(^)e-'' (2) 

with = 0.0035(1 + z)-2-20^ a = j^glog(Mg,i)-M,^ ^ 

Il.16-f0.17z-0.07z2, anda(z) = -1.18-0.082z. A galaxy 
mass, Afgai, in units of M©, is drawn from this distribution 
(Eq. ^ in the range 7 < log(Mgai) < 12. For galaxies at 
redshifts beyond z = 4 we use the above distribution eval- 
uated at z = 4, i.e., we assume no evolution in the galaxy 
mass distribution at earlier times. This assumption re- 
flects the lack of information on galaxy mass distributions 
at high redshift. Although this may be a decent approxi- 
mation at intermediate redshifts (4 ^ z ^ 10), it is almost 
certainly incorrect at very high redshift, where only low 
mass galaxies (and Population HI stars) are forming. 

Galaxy Metallicity. We assume that the average oxy- 
gen to hydrogen number ratio (Fon) of a typical galaxy 
depends on its mass as: 

log(i^oH) = log(10i2O/H) = 

sz + 1.847 log(Mgai) - 0.08026(log(A/gai))2 

with a redshift independent normalization sz — —1.492 
(Tremonti et al. 2004). It has been suggested that the 
functional form of this mass-metallicity relation is redshift 
independent (Erb et al. 2006; Young & Fryer 2007), with 
only the normalization factor, sz, varying with redshift. 
We characterize the redshift dependence of galaxy metal- 
licities using the average metallicity relation from Pei, Fall, 
& Hauser (1999): 

( IQ-"^^ z < 3.2 
Z oc \ io-bi-b2z 3.2 < z < 5 , (4) 
[ IQ-'^i-^^^ z > 5 

which implies evolution of the normalizing factor, sz, with 
redshift given by 

C -a2Z - 1.492 z < 3.2 

sz oc < ~b2Z - 3.2(a2 - 62) - 1-492 3.2 < z < 5 

[ -C2Z - 5(62 - c2) - 3.2(a2 - 62) - 1.492 z > 5 

With this approach we assume that the oxygen abundance 
(used in -Fqh) correlates linearly with the average abun- 
dance of elements heavier than Helium (which is what is 
provided in the metallicity measure, Z). The normal- 
ization given by Pei et al. (1999) is 02 — 0.5, bi = 
0.8, 62 = 0.25, ci = 0.2, C2 = 0.4. We also derive 
the coefficients from the metallicity dependence proposed 
by Young & Fryer (2007), based on uhraviolet (GALEX, 
Sloan Digital Sky Survey), infrared (Spitzer), and neu- 
trino (Super-Kamiokande) observations (Hopkins & Bea- 
com 2006): 02 = 0.12, bi = -0.704, 62 = 0.34, a = 
0.0, C2 = 0.1992. The Pei et al. (1999) normalization re- 
sults in a relatively slow metallicity evolution as compared 
with the model of Young & Fryer (2007), which allows for 
a rapid increase of the average metallicity of stars with 
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time. This is shown in Figure [T] (top panel) , which pro- 
vides the metaUicity distribution of stars. It is apparent 
that for the fast model the metaUicity of stars is on aver- 
age an order of magnitude higher than for the slow model 
at the same high {z = 7) redshift. 

Using Equations[3]and[5l we estimate the average galaxy 
oxygen to hydrogen ratio, given the mass and redshift of 
a galaxy. We express our estimators relative to solar: 

^OH,gal ^on/i^OH,©, (6) 

where Foh,© = 4.9 x 10^ is the solar value (Tremonti et 
al. 2004). ' 

Galaxy Stellar Population. We use the above galaxy 
metaUicity estimator to formally delineate different stellar 
populations: 

^oH,gal < 10"^ Population III 

10-^ < ^'oH,gal ^ 10-1 Population II . (7) 

^OH gal 10^^ Population I 

Our choice for the lower metaUicity bound (10^'*) on Pop- 
ulation II marks the point where metals are abundant 
enough to provide sufficient cooling in the collapse of gas 
clouds, and thus star formation significantly deviates from 
the Population III stage (e.g., Mackey, Bromm, & Hern- 
quist 2003; Smith et al. 2008). It also marks the point 
where winds of massive stars are sufficiently strong to pre- 
vent the formation of pair-instability supernovae (Heger et 
al. 2003). The choice for the Population II upper bound 
is dictated by the observation of Population I and II stars 
in the Milky Way (e.g., Binney & Merrifield 1998; Beers & 
Christlieb 2005). Population I stars (e.g., disk and bulge 
stars) have metallicities that are approximately solar, with 
variations of a factor of five. Population II stars (e.g., halo 
populations, such as stars in globular clusters) have metal- 
licities at or below lO-^-lO"^ solar. 

Each galaxy is assumed to host just one stellar popula- 
tion, determined by its metaUicity. Since we draw a large 
number of galaxies at any given redshift (see Eq. [2]) via 
Monte Carlo simulations, the use of an average metaUicity 
measure is appropriate. 

In Figure [T] (bottom panel) we show our assumed star 
formation history as a function of redshift, with the con- 
tribution from Population II stars highlighted. Popula- 
tion II stars form over a wide range of redshifts, starting 
at z ~ 22, and peaking at z ~ 7 for the fast metaUicity 
evolution model. For the slow metaUicity evolution, the 
interstellar medium is enhanced with metals much later, 
and Population II stars begin forming only at redshifts 
z ~ 11, peaking at z ~ 3. Although the metaUicity evo- 
lution is poorly constrained, the two adopted models are 
most likely extremes, with the true evolution somewhere 
in between. For example, the transition from Population 
III to Population II stars is found to occur at z ^ 11 and 
z 22 for our slow and fast metaUicity evolution models, 
respectively. More sophisticated models place this limit 
at redshifts z - 15-20 (e.g., Mackey et al. 2003). Re- 
cent results show that Population III stars can form as 
recently as z ~ 5 (Schneider et al. 2006; Tornatore, Fer- 
rara, Schneider 2007), but with rates significantly below 
that of Population II stars. 



Population III stars were forming prior to the measured 
redshifts of GRB 080913 and 090423. Individual Popu- 
lation III stars may give rise to long GRBs in coUapsar 
models with enhanced rotation (e.g., Yoon, Langer, & 
Norman 2006). However, as the lifetime of massive Pop- 
ulation III stars is very short (~ 3-6 Myr; e.g., Schaerer 
2002), these GRBs would be observed at much higher red- 
shifts (z ^ 10), tracing the Population III star forma- 
tion rate with little delay. If binary Population III stars 
are formed, they could evolve into NS-NS/BH-NS systems 
that in turn, after an appropriate delay before merging, 
would give rise to short GRBs at potentially much lower 
redshift (although see Ripamonti & Abel 2004 and O'Shea 
& Norman 2007, both of whom find that fragmentation 
does not occur in pristine clouds, implying that no Pop- 
ulation III binaries would form). It has been argued that 
GRBs from such binaries, if they do form, would result in 
negligible detection rates for current detectors, including 
Swift (Belczynski et al. 2007). 

The transition from Population II to Population I star 
formation begins at z ~ 7 for the fast or at z 3 for 
the slow metaUicity evolution model, continuing gradually 
until the present. Regardless of the metaUicity evolution 
model. Population I stars are the dominant contribution 
to the star formation rate below z ~ 2, although a small 
number of Population II stars are also predicted at these 
lower redshifts. The earliest Population I stars begin form- 
ing at z '-^ 7, but this is a negligible fraction of the total 
star formation at this redshift, even for the extreme fast 
metaUicity evolution model (see Fig. [TJ top panel). Signif- 
icant formation of Population I stars is found for redshifts 
z <J 6 (fast) and z <^ 2.5 (slow), and thus Population I 
stars are unlikely to be the progenitors of high-z bursts. 

Regardless of metaUicity evolution models. Population 
II stars dominate the star formation rate at redshift 6 <^ 
z ^ 10 (see Fig. [1]), and are the likely progenitors for both 
GRB 080913 and GRB 090423. Therefore, in the following 
we focus exclusively on Population II stars. 

3. STELLAR EVOLUTION: MODEL/RESULTS 

Evolutionary Code. We use the StarTrack population 
synthesis code (Belczynski et al. 2002, 2008a), which em- 
ploys the stellar evolution models of Hurley, Pols, & Tout 
(2000), to calculate the population of low metaUicity stars, 
with Z = 0.0001 (- 10^2.32-^)^ gyoi^g j^^.^^ = 2 X 10^ 

massive binary systems, with the mass of the primary (the 
initially more massive) member. Mi, drawn from a power- 
law initial mass function (IMF) with exponent —2.7 in the 
range 5-150 M©. The secondary member's mass, M2, is 
drawn randomly from a uniform initial mass ratio distri- 
bution {q — M2/M1 £ — 1), where we only evolve binaries 
that have secondaries with M2 > 3 Mq . Initial orbital sep- 
arations, a, are drawn from a uniform random distribution 
in In a (i.e., pdf oc 1/a; Abt 1983), and eccentricities are 
taken from the distribution pdf = 2e (e.g., Heggie 1975; 
Duquennoy & Mayor 1991). 

We also evolve a population of iVgin = 2 x 10^ massive 
single stars, with initial masses chosen as described above 
for the primary members of the binaries. The lower mass 
limits are chosen such that all stars that can potentially 
form either neutron stars or black holes are included in the 
calculations. 
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If the mass range is extended to the hydrogen-burning 
hmit (~ 0.08 M0), one can estimate the total stellar mass 
that corresponds to our calculation (i.e., we evolve only 
the fraction of mass that is formed in massive stars, but 
for the purpose of rate calculations we also estimate the 
total stellar mass). We use a 3-component power law IMF; 
for stars with masses 0.08-0.5 the exponent is —1.3, 
for stars in the mass range 0.5-1 M© we use —2.2, and for 
masses above 1 Mq we use —2.7 (Kroupa, Tout, & Gilmore 
1993; Kroupa & Weidner 2003). Note that our population 
assumes an equal number of single and binary star sys- 
tems, corresponding to a binary fraction /bi = 50% (or a 
binary system membership probability of 2/3 for any ran- 
domly selected star). Spectroscopic studies demonstrate 
that binary fractions are significant, especially for mas- 
sive stellar populations (e.g., > 70% for the Westerlund- 
1 cluster (Clark et al. 2008)7 or ^ 50% for NGC6231 
(Sana et al. 2008)). The total mass in our simulation, 
in single and binary stars over the entire mass range, is 
Msin, = 6.1 X 10^ Mq. 

Short GRBs. We extract from the population synthesis 
sample the NS-NS and BH-NS binaries with delay times 
shorter than 15 Gyr. The total number of short GRB pro- 
genitors is found to be rtshort — 7.6 x lO'^ with 4.6 x 10'^ 
NS-NS and 3.0 x 10^ BH-NS systems. The delay time, 
idcb is the time a given binary takes from the zero age 
main sequence to form a double compact object (evolu- 
tionary time, tevoi) plus the time a double compact object 
takes to merge due to emission of gravitational radiation 
(merger time, ^merger)- The resulting delay time distribu- 
tion is shown in Figure [2l The median delay time is 0.1 
Gyr, with a mean of 1.5 Gyr and a standard deviation of 
2.9 Gyr. Contrary to conventional wisdom, most mergers 
are expected shortly after the stars form, with only a small 
fraction (3.7%) of delay times in excess of 10 Gyr (qualita- 
tively similar results for Population I stars were presented 
in Belczynski et al. 2006). 

Long GRBs. We extract black holes (formed in both 
isolated stars and in binaries) that (i) formed through di- 
rect collapse (e.g., Fryer & Kalogera 2001) and (ii) formed 
from a progenitor that at the time of core collapse was 
hydrogen depleted. The first condition ensures rapid ac- 
cretion (~ 0.1 — IMq s~^) onto a black hole. Popham et 
al. (1999) found that at least some of the power sources 
for GRBs are much more efficient at higher accretion rates, 
and Young & Fryer (2007) argue that current observations 
are better fit by these direct-collapse black holes. For the 
adopted Population II metallicity, single stars with initial 
mass Mzams > 23.8 Mq form black holes through direct 
collapse. Note that the limiting mass was 40.0 for so- 
lar metallicity, as obtained in Fryer & Kalogera (2001). 
Here we use different evolutionary models (Hurley et al. 
2000) and we employ stellar winds that scale with metal- 
licity as oc Z^^^, and find that stars with lower mass can 
potentially form black holes directly. However, we note 
that had we used 40.0 M0 as a limiting mass for direct 
black hole formation, our predicted rates for long GRBs 
would have remained virtually the same since condition 
(ii) selects stars with mass that are above or close to the 
previous limit (see below). Our limit, which takes into 
account the effect of metallicity and stellar winds, is fur- 
ther modified by the presence of a binary companion (from 



rejuvenation and/or mass loss in close binaries). 

The second condition allows for a GRB jet to punch 
through the outer layers of a star. A star may lose its 
H-rich envelope either through binary interactions or via 
stellar winds. Given the employed single star models, and 
the adopted metallicity, stars with Mzams > 35.9 lose 
their entire H-rich envelopes via stellar winds. Binary in- 
teractions, on the other hand, may remove the envelope 
of a star of arbitrary mass. These conditions are insuffi- 
cient to produce a GRB, since there must be just enough 
angular momentum in a collapsing star to form a lasting 
accretion torus (e.g., MacFadyen & Woosley 1999; Podsi- 
adlowski et al. 2004) . As angular momentum transport is 
not taken into account in this study, our predictions only 
indicate potential long GRB progenitors, i.e. most likely 
only a (small) fraction of the events that satisfy (i) and 
(ii) will actually produce long GRBs. 

The total number of long GRB progenitors is found to 
be niong = 225 x 10^ with 103 x 10^ and 122 x 10^ di- 
rect black holes (DBHs) from H-depleted progenitors from 
single and binary stars, respectively. The delay time dis- 
tribution is displayed in Figure [21 This represents the 
evolutionary time (iovoi) of a given star to form a black 
hole. The delay times for long GRB progenitors are very 
short, ^ 3.5-12 Myr (as black holes form from massive and 
short-lived stars), with the median of the distribution at 
5.1 Myr, a mean at 5.3 Myr, and a standard deviation of 
1.1 Myr. These events thus trail the star formation rate, 
SFR(z), with very little delay, so that essentially one can 
assume that Rate(DBH) oc SFR until very large redshifts. 

4. COSMOLOGY: MODEL 

We adopt a flat cosmology with Hq = 
70.0 km s"^Mpc"\ = 0.3, and Ha = 0.7 (and thus 
Ok = 0). Small (^ 10%) changes in the cosmological 
parameters leave our results essentially unaltered. 

For clarity in what follows, we briefly review the rele- 
vant cosmological expressions. The relationship between 
redshift and time is given by 



t{z) = in 



dz' 



{1 + z')E{z')'^ 



(8) 



where in — ^/Hq — 14 Gyr is the Hubble time (e.g., Hogg 
2000), and E{z) = ./17m(1 + + f7k(l + z)^ + flA). 
The resulting age of the Universe is t{0) — 13.47 Gyr. 
It is to be noted that this is the restframe time, consti- 
tuting what is measured by the wristwatches of local ob- 
servers at any time and place in the Universe, and thus 
is the appropriate quantity to use when discussing local 
physical processes. Observed event rates from sources dis- 
tributed over cosmic time must be corrected for time dila- 
tion, R{0) = R{z)/{l + z). 

The comoving volume element, dV/ dz, for a given solid 
angle, dil, is given by 



dz Ho E{z) 



dn. 



(9) 



where c is the speed of light in vacuum, and where the 
comoving distance, D^, is given by 



-no Jo 



dz' 

E{z' 



(10) 
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There are a number of steps necessary to calculate the 
expected GRB event rates. To summarize: for each bin 
in time, we estimate the total star formation rate. We 
then calculate the fraction of these stars which are Pop- 
ulation II by incorporating the metallicity distribution of 
the galaxies at the given time. We then employ population 
synthesis results to determine the fraction of Population II 
stars that are progenitors of GRBs. With these results, we 
are able to apply the appropriate delay times and trans- 
form the progenitor formation rate in each bin in time to 
an actual GRB event rate, distributed over later times. 
We repeat this procedure over all bins in time, which re- 
sults in a total GRB event rate as a function of time. We 
now go through each of these steps in detail. 

We bin time over the range 0.13 < t < 13.47 Gyr 
{0 < z < 25), with a constant width of dt = 10 Myr. For 
each bin we obtain the star formation rate using Equa- 
tion [U evaluated at the center of the bin. We then gen- 
erate a Monte Carlo sample of 10"* galaxies in the bin, 
distributed according to Equation [21 with a total mass in 
galaxies given by Afgai.tot- For each galaxy we estimate 
its average metallicity from Equation [3] and, based on our 
population criterion (Eq. [7|), we determine the number of 
galaxies containing Population II stars. The total mass of 
galaxies with Population II stars is denoted Afgai.n. The 
fraction of mass that has formed Population II stars in a 
given bin is thus given by 



sal, II 



pop, II 



gal, tot 



(11) 



since we have assumed that each galaxy hosts only one 
stellar population (determined by its average metallicity; 
see §2). 

We measure rate densities in units of yr~^ Gpc"'^, in 
comoving units (time and space). In our scheme a delay 
for a given event constitutes a shift in time of the relevant 
fraction of the rate density (see Eq. [T^ . The population 
synthesis results yield a list of Population II GRB progeni- 
tors, and their corresponding delay times, and we use all of 
them in each time bin. Thus, for each bin we have a total 
of rishort = 7.6 X 10'^ and niong = 225 x 10^ progenitors for 
short and long GRBs, respectively. Each progenitor in our 
list now represents a fraction of the total star formation 
rate density in the given bin that generates GRBs. This 
fraction is given by: 



/si 



pop, II 



A/si 



sfr(t), 



(12) 



where Mgim is the total stellar mass corresponding to our 
population synthesis simulations (see § 3) , and fsim has 
units of yr~^ Gpc^'^. 

For each progenitor in the list we choose a random start- 
ing time, ^0, within the given time bin, and propagate the 
progenitor formation rate density, /simj forward in time to 



incw — ^0 ^dcl. 



(13) 



The delay time, tdcii marks the time elapsed from the for- 
mation of the star to the potential GRB event, as described 
in § 3, and t^cw is the time at which a given GRB event ac- 
tually occurs. We thus convert the formation rate density 



of progenitors into GRB event rate densities. By repeating 
the above series of steps for each time bin, we arrive at a 
total comoving GRB event rate density, n,:cst{t)- 

5. INTRINSIC RATES 

In this section we calculate the intrinsic rate of short 
and long GRBs from Population II progenitors. In the 
first subsection we discuss correction factors due to model 
uncertainties, while in the second subsection we provide 
estimates of the intrinsic GRB rates. 

5.1. Correction factors: population synthesis 

Both the direct black hole formation and double com- 
pact object merger rates obtained from population synthe- 
sis are only first-order approximations for long and short 
duration GRBs, respectively. For example, the only con- 
straints placed on our long-duration GRB progenitors were 
the assumption that the black hole formed from direct col- 
lapse and that the star, at collapse, was a He star. Obser- 
vations of nearby GRB/SN associations suggest that the 
progenitor must lose most of its Helium, because, to date, 
every supernova observed associated with a GRB is a type 
Ic supernova (see Fryer et al. 2007 for a recent review). We 
have not placed additional constraints requiring the loss 
of the He envelope. In addition, we have not placed con- 
straints on the angular momentum profile, which is likely 
to play an important role in GRB formation. 

The type Ic constraint poses a problem for many of the 
current progenitor scenarios and, hence, makes it difficult 
to predict a reliable rate from population synthesis studies 
alone. Ignoring this constraint, we can make only rudimen- 
tary estimates of the GRB event rates. With the rotation 
rates suggested by Yoon et al. (2006), roughly 10% of the 
collapsing stars in our simulations form GRBs. Given our 
rate estimates for both long and short bursts, this indi- 
cates that the odds are roughly 10 to 1 that GRB 080913 
and 090423 were long-duration bursts. However, Podsi- 
adlowski et al. (2004) have argued that only 1% of all 
black hole systems would form GRBs. Although these au- 
thors mention that this fraction could potentially increase 
at low metallicity, if we keep their 1% figure we find a 
long-duration GRB rate 100 times smaller than the one 
we calculated for the direct black hole formation channel. 
Thus, the likely rate for long-duration GRBs is reduced by 
a factor of esyn,iong ^ 0.01-0.1. 

Our rate for short GRBs is also an over-estimate. We 
include all BH-NS mergers, while neglecting system config- 
uration (inclination/mass ratio) and black hole spin (only 
1-40% percent of mergers are expected to form BH torus 
configurations and thus potentially result in a GRB: Bel- 
czynski et al. 2008b). If GRBs indeed require the forma- 
tion of a black hole (e.g., Janka & Ruffert 1996; Oechslin 
& Janka 2006; Lee & Ramirez- Ruiz 2007), the fraction of 
NS-NS mergers is limited by the maximum neutron star 
mass. Most known NS-NS mergers will form central com- 
pact objects with masses Mc <^ 2.5 M0 (Belczynski et al. 
2008c). The maximum neutron star mass must be less 
than this value, or these mergers will not form black hole 
accretion disks. However, the requirement of a black hole 
is a weak constraint as neutron star accretion disks as well 
as proto-magnetars may also produce short GRBs (Usov 
1992: Kluzniak & Ruderman 1998; Dai et al. 2006; Met- 
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zger, Quataert & Thompson 2007). The low metaUicity 
population of double compact objects consists of ~ 60% 
of NS-NS binaries and - 40% of BH-NS systems (see § 3). 
Thus the reduction due to the required configuration of 
BH-NS systems does not lead to a significant decrease of 
short GRB rates. The likely rate for short-duration GRBs 
is thus reduced by approximately Cgyn. short ^ 

0.6-^0.8. 

Additionally, the formation rates of double compact ob- 
jects as calculated in population synthesis are subject to 
uncertainties that are associated with some poorly under- 
stood evolutionary processes that are important in the 
evolution of massive binaries (e.g., winds, mass transfer 
episodes, natal kicks). Belczynski et al. (2002) presented 
a comprehensive parameter study of double compact ob- 
ject formation rates. If the most unrealistic and unphys- 
ical models are omitted, it was shown that NS-NS and 
BH-NS merger rates vary at most by factor of ~ 6 (up 
and down) from the reference model (best guess for evolu- 
tionary parameters) value. We use this as an estimate of 
the population synthesis uncertainty in the rates and thus 
renormalize correction factors to Csyn, short ~ 0.1-4.8. 

In summary, the population synthesis or model related 
correction factors are 

Csynjong -0.01-0.1 (14) 
Csyn.short = 0.1-4.8. (15) 

5.2. Intrinsic rate estimate 

Our estimate of the intrinsic GRB event rate density in 
the rest frame is given by 

A/'rcst(t) = Esyn J^rcst (0 Y^'^ Gpc"^, (16) 

where Csyn was introduced in the previous subsection, and 
rirest was introduced in § 4. The rate density of GRBs 
that would be observed is appropriately time dilated: 
A/'obs(i) = A/'rest(i)/(l + z{t)). In a slightly confusing ac- 
cident of notation, the units of both of these quantities is 
yr~^ Gpc~'^, where in one case years are measured in the 
rest frame of the GRB population, and in the other case 
years are measured by a clock at the local observatory. 

It is to be noted that some events have such long de- 
lay times, or are associated with such recent star for- 
mation, that the GRB events will happen in our future 
(inew > 13.47 Gyr), and thus are not included in the event 
rate density estimates. We note that this affects only short 
GRB progenitors, as delay times for long GRB progenitors 
are very short (see Fig. [2]). 

We are also interested in the intrinsic GRB event rate, 
in addition to the event rate density. We present intrinsic 
event rates in the observer frame. To determine the rate 
to some specified redshift, z, we integrate the GRB event 
rate density (modified by time dilation) over the comoving 
volume: 

where we have integrated over the entire sky (J dil = Att). 
It is important to note that, although the above integral 
is in redshift, it is integrating out the volume (given by 
Eq.[5]), and not the time: the original units of rate density 



(yr ^ Gpc ^) have been converted into a rate (yr ^) and 
not into a density (Gpc~'^). A^(< z) represents the rate 
in the observer frame of all GRB events (bursts per year) 
integrating to a limiting redshift of z. Since N denotes the 
intrinsic rate, we have not yet applied observational selec- 
tion effects due to such things as beaming, brightness, or 
instrumental response (these are calculated in §6.1). 

The population synthesis code naturally distinguishes 
between long and short GRBs, and we carry this identi- 
fication through our calculations, quoting results for both 
long (Mong) and short (A^short) GRBs. The resulting in- 
trinsic GRB rate densities and event rates are shown in 
Figures [31 [H and [S] As expected for long GRBs, these 
rates closely track the star formation history of Popula- 
tion II stars: the majority (90%) of events are found at 
redshifts z ^ 3-15 with a peak at z ~ 7 (fast metallicity 
evolution model), or at z ~ 2-11 with a peak at z ~ 3 
(slow model). Short Population II GRBs start appearing 
slightly later than long Population II GRBs, and their dis- 
tribution extends to lower redshift (z ~ 0). However, the 
peak of the short GRB rate is near the long burst peak, 
at z ~ 7 (fast) and z ~ 3 (slow), due to the large abun- 
dance of very short delays. We note that even a delay as 
short as 0.1 Gyr is sufficient to result in a many kpc dis- 
tance between the site of star formation in the host galaxy 
and the eventual merger location in the halo (e.g.. Fryer, 
Woosley, & Hartmann 1999; Bloom, Sigurdsson, & Pols 
1999; Belczynski et al. 2006). The short GRB progeni- 
tors considered here (NS-NS/BH-NS mergers) have a tail 
of very long delay times (median 0.1 Gyr, and standard 
deviation 2.9 Gyr), resulting in substantial intrinsic rates 
of Population II mergers at low redshift. 

We note that in the redshift range 6 < z < 8, the 
intrinsic rates for long and short GRBs are comparable 
(see Fig. [S]). At redshift z = 6.7, the rates are Mong = 
10 X 10^-153 xlO^ yr-i and A^short = 3 x 10^-141 x 10^ yr^^ 
per unit redshift, in the observer frame. At redshift z = 8.1 
these rates become Mong = 9 x 10^-91 x 10^ yr^^ and 
A'^short = 2 X 10^-84 X 10'^ yr~^ (the ranges encompass both 
of our metallicity evolution models). The above similarity 
of rates arises for two reasons: first, the population syn- 
thesis rates of short (esyn, short 'T'short = 0.8 X 10'^-36.5 X 10'^) 
and long (esyn.iongfiiong = 2.3 x 10^-22.5 x 10^) GRB pro- 
genitors are comparable; second. Population II star for- 
mation rates for the two metallicity evolution models are 
almost the same (~ 0.1 Mq yr~^ Mpc~'^; see Fig. [T|) at 
these redshifts. 

For comparison, we have also calculated the rate of su- 
pernovae utilizing the same procedure presented above. 
The cumulative rate of core collapse supernovae (Type 
II and Ib/c) is estimated to be ~ 8 s~^ (integrated to 
redshift z = 10), which is comparable to the rates from 
recent empirical estimates of SFR(z), e.g. Hopkins and 
Beacom (2006). 



6. OBSERVED RATES 

In this section we predict the Swift detection rate of 
short and long GRBs from Population II progenitors. In 
the first subsection we discuss observational selection ef- 
fects, while in the second subsection we provide estimates 
of the observed GRB rates. 
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6.1. Correction factors: observational selection effects 

Not all high-redshift GRBs will necessarily be observ- 
able. There are many factors which impact whether or 
not a given GRB will be observed, including the lightcurve 
and spectral energy distribution of the burst (both of them 
redshifted), as well as the spectral and time sensitivity of 
the detector. A precise calculation of the selection func- 
tion of GRBs is well beyond the scope of this paper. The 
GRB population is tremendously heterogeneous, and the 
triggering algorithms of an instrument such as the BAT 
on Swift are highly complex (Band 2003, 2006; Fenimore 
et al. 2004). In what follows we make a number of sim- 
plifications, and derive an approximation for the observed 
rate of short and long GRBs. 

We assume that all GRBs have the same distribution of 
average isotropic luminosity, Liso, given by a Gaussian in 
log(Liso), centered on Ljso = 10^° erg s~^, with Itr = 1 dex. 
We assume all short GRBs have a duration of tgo = 0.3 s, 
while all long GRBs last tgo = 30 s. The distributions of 
total energy (integrated flux; calculated as £ = Lisoigo) 
for short and long bursts are presented in Figure [SI 

We approximate the Swift GRB trigger (for the BAT 
instrument) as a_fluenccE] threshold. We take the long 
GRB trigger as £iong = 3 x 10^" erg for a long burst at 
z = 1 (see, e.g., Fig. 2 of Kistler et al. 2008 and Fig. 2 
of Butler et al. 2007). In other words, a long GRB at 
z = 1 must emit a total isotropic-equivalent energy satis- 
fying £" > 3 X 10^° erg to be observable. The corresponding 
short GRB threshold is given by ^short = 10^" erg at z — 1. 
The difference between the thresholds is chosen to mimic 
some of the complexity in the BAT triggers, including in- 
tegration and readout times, as well as sensitivity to the 
differing spectral energy distributions (SEDs) of long and 
short bursts (Fenimore et al. 2004; Band 2003, 2006). For 
a burst at redshift z, we draw an energy from the relevant 
distribution shown in Figure[6l and then calculate whether 
the fluence of the burst satisfies our trigger threshold cri- 
terion: 

£>{DL{z)/DLiz^l)f£, (18) 



where Dl{z) is the luminosity distance out to redshift z. 
Only those bursts which satisfy this inequality are taken 
to be observed. The resulting fraction of long and short 
bursts observed with Swift, esed, as a function of redshift, 
is plotted in Figure [6l 

Another factor we have not yet considered is beaming, 
and the resulting visibility of GRBs. In general short GRB 
jets appear to be wider, with opening angles of 9 > 10° 
(e.g.. Fox et al. 2005; Berger 2007; Metzger, Piro, & 
Quataert 2008) as compared with long GRBs, 0^5° 
(e.g.. Frail et al. 2001; Bloom, Frail, & Kulkarni 2003; 
Grupe et al. 2006; Soderberg et al. 2006). This would 
imply that long GRBs are ^ 4 times less likely to be de- 
tected, leading to further reduction of the observed rates 
by factors of 

Cbcam.long = 0.002 (19) 



short 



0.008 



(20) 



for long and short GRBs, respectively. 



6.2. Observed rate estimate 

We now provide a prediction for the Swift detection rate 
for Population II GRBs. We express rates per unit redshift 
to emphasize the redshift range relevant for GRB 080913 
and GRB 090423. The predicted detection rate is given 
by: 

= /Swift eSED Cbcam 



dz 



dz 



yr 



(21) 



where N is the intrinsic rate of GRBs (see § 5.2), ebcam is a 
correction factor for GRB beaming (see §6.1), esED is the 
fraction of GRBs above the Swift detection threshold (see 
Fig. Eland §6.1), and /swift = 1.4/(47r) = 0.1 represents 
the Swift sky coverage (Gehrels et al. 2004; Barthelmy et 
al. 2005). We plot our prediction for the observed rate of 
both short and long GRBs in Swift in Figure [T] In the red- 
shift range z = 6.2-7.2 centered on z = 6.7 (GRB 080913), 
we expect Swift GRB detection rates of: 



ift,long 



Swift, short 

For the redshift 
(GRB 090423), 
come: 

^ Swift., \o\\g 

^ Siui/t, short 



0.33-5.0 yr"' 
0.01-0.70 yr' 



at z = 6.2-7.2 
^ at z = 6.2-7.2. 



(22) 
(23) 



range z = 7.6-8.6 centered on z = 8.1 
the predicted Swift detection rates be- 



0.21-2.2 yr^i at z = 7.6-8.6 (24) 
0.004-0.24 yr^i at z = 7.6-8.6. (25) 



These results suggest that GRB 080913 and GRB 
090423 are most likely associated with the deaths of mas- 
sive stars (long GRBs), rather than double compact ob- 
ject mergers (short GRBs). At redshifts z = 6.7 (GRB 
080913) and z = 8.1 (GRB 090423) long GRBs are ex- 
pected to be ~ 10 times more frequent than short GRBs 
in the Swift sample. We emphasize that we have taken an 
optimistic value for the short burst observational thresh- 
old; the true short GRB rate may be up to an order of 
magnitude smaller (see §6.1 for details). Our results thus 
represent an upper-limit to the short GRB fraction at these 
high redshifts. Although the intrinsic event rates for col- 
lapsars and double compact object mergers are similar at 
these high redshifts, the observational difference is due to 
the very distinct observational selection effects between 
these two event types. We note that our results are consis- 
tent with Zhang et al. (2009), who classified both bursts as 
coUapsars based upon their observational properties. We 
have ignored recent evidence that the long GRB popula- 
tion evolves with redshift, increasing the number of ob- 
served bursts at high redshift. In this sense, our results 
can be considered a lower limit to the long GRB rate (al- 
though we note that the corresponding evolution of short 
bursts is currently poorly constrained). Further details of 
the nature of GRB intrinsic properties and observational 
selection effects can be found for example in Guetta, Pi- 
ran & Waxman 2005; Guetta & Piran 2006; Berger 2007; 
Kistler et al. 2008; Salvaterra et al. 2009 and Virgih, 
Liang & Zhang 2009. 



^Our fluence threshold can easily be rephrased as a trigger on peak flux/luminosity, with all of the units being converted from [erg] to [erg 
s~^]. The results are identical. In addition, we are calculating the fluence from an isotropic luminosity, but the results are unchanged if this is 
taken to be a beamed luminosity. 
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7. SUMMARY 

We have calculated the evolution of stars at high red- 
shift, exploring possible progenitors for GRB 080913 (at 
z = 6.7) and GRB 090423 (at z = 8.1). We find that in the 
redshift range 6 ^ z <^ 10 the majority of stars are Pop- 
ulation II (low metallicity). By this time Population III 
stars (metal free) have already finished their evolution, and 
Population I stars (metal rich) arc only just beginning to 
form. We have adopted two extreme models for the metal- 
licity evolution history, which most likely bracket the true 
history; our final results are insensitive to the choice of 
model. We find that the progenitors of GRB 080913 and 
GRB 090423 are likely to have been Population II stars. 

For Population II stars we subsequently calculated the 
rates of double compact object mergers (NS-NS and BH- 
NS) and the rates of coUapsars (deaths of massive stars). 
The former are thought to be associated with short GRBs, 
while the latter are believed to be responsible for long 
GRBs. We find that the intrinsic rates for short GRBs 
(mergers) and long GRBs (coUapsars) are comparable at 
high redshift (see Fig. [5]) . However, the observational se- 
lection effects (beaming, differing intrinsic fiuence, and in- 
strumental response) make long bursts more likely to be 
seen by a satellite like Swift. Our main result is shown 



in Figure [3 where we plot our predicted Swift detection 
rates for both short and long GRBs. On average, the de- 
tection rates for long GRBs are 10 times higher than the 
rates for short GRBs. At the redshifts of GRB 080913 and 
GRB 090423, the rates are - 1 yr^^ and - 0.1 yr~i per 
unit redshift, for long and short GRBs, respectively. We 
emphasize that our calculations only consider Population 
II stars, as these stars dominate the star formation rate at 
redshift 6 <^ 2; ^ 10 (see Fig. [Ij, and are thus appropri- 
ate for the two highest known redshift GRBs. In future 
work we will include all generations of stars. 

We conclude that both GRB 080913 and GRB 090423 
are most likely to have been long bursts, resulting from 
the deaths of massive stars (coUapsars) from Population 
II progenitors. With the observation of GRB 080913 and 
090423, gamma-ray bursts have entered a unique high- 
redshift regime. As future data fills out the high- 2; tail of 
GRBs, these systems will become one of the most powerful 
probes of the star formation, stellar death, and chemical 
enrichment of our Universe. 

We thank the anonymous referee for encouraging us to 
incorporate observational selection effects into our analy- 
sis. KB acknowledges support from KBN grant N N203 
302835. 
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Fig. 1. — Bottom panel: Star formation rate history. We have adopted an extinction corrected SFR model from Strolger 
et al. (2004): all stars. We also show the rate we have obtained for just Population II stars (see §2): the population is 
rather different for slow (solid) and fast (dashed line) metallicity evolution model. Top panel: Metallicity distribution of 
stars at redshift z — 7.0. The results are shown for the two different adopted metallicity evolution histories: fast (Young & 
Fryer 2007) and slow (Pei et al. 1999). We mark our definition for Population II stars: it is clearly seen that independent 
of metallicity evolution majority of the stars belong to Population II at redshift that is typical for stars that may become 
(after appropriate delay) progenitors of GRB 080913. 
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Fig. 2. — Top panel: Delay times for long GRB progenitors: direct black holes that are formed out of H-depleted stars 
(both single and binary). Note the very short delay times of <; 6 Myr. Delay is evolutionary time a star takes from the 
formation to a collapse. Bottom panel: Delay times for short GRB progenitors: NS-NS and BH-NS mergers. Note that 
these events have significantly longer delay times than for the long GRB progenitors, with a median of 0.1 Gyr. Delay 
time is evolutionary time plus merger time. 
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Fig. 3. — Short and long GRB intrinsic event rate densities (in the restframe) originating from Population II stars for 
slow and fast metallicity evolution. These are expressed per year, per cubic comoving Gpc. 
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Fig. 4. — Intrinsic event rate of short and long GRBs that originate from Population II stars, per year (in the observer 
frame), as a function of the depth of the survey in redshift, z. The rates are shown for slow and fast metallicity evolution 
models. 
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Fig. 5. — Intrinsic event rate of short and long GRBs that originate from Population II stars, per year (in the observer 
frame) presented per unit redshift. These plots show the derivative of curves presented in Figure [H highlighting the 
redshift range in which GRB events from Population II stars contribute. 
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Fig. 6. — Top panel: Distribution of total energy {£ — L{soHo) for short and long GRBs. Short and long GRBs are 
assumed to have the same distribution of average isotropic luminosity, Liso, described by a Gaussian in log(Liso) with 
mean at log(_Liso/[erg s~^]) = 50 and a = 1.0. A single typical duration time was assumed for all short (igo = 0.3s) 
and long (igo = 30s) GRBs. The adopted Swift detection thresholds (at z = 1.0) are also shown; £ — 10^°, 3 x 10^°erg 
for short and long GRBs, respectively. See the text for details. Bottom panel: The fraction of GRBs above the Swift 
detection limit, esED, as a function of redshift. Note that at redshift z — 6.7 (GRB 080913) the fraction of detectable 
long GRBs is esED,iong — 0.162, while for short events it is only csed, short = 0.006. At redshift z = 8.1 (GRB 090423) the 
difference in detection probability is even larger: esEDjong = 0.118 and csed, short = 0.004. 
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Fig. 7. — Predicted Swift detection rates of short (NS-NS/BH-NS mergers) and long (collapsars) GRBs. Rates are for 
GRBs originating exclusively from Population II stars, and hence can only be compared to Swift data in the redshift 
range in which Population II stars are the dominant stellar population: 6 z 10 (as marked on the plot). For low 
redshifts (z <^ 6) there is a significant contribution to GRB rates from Population I stars (not shown here), and for very 
high redshifts {z ^ 10) there may be a contribution from GRBs originating from Population III stars (not shown here). 
Note that GRB 080913 and GRB 090423 are both found in the redshift range in which the majority of GRB events are 
predicted to originate from Population II stars. Furthermore, our results favor the coUapsar origin for these GRBs. The 
two panels show the rates for a fast (early metal mixing into stars; top panel) and a slow (late mixing; bottom panel) 
metallicity evolution model. In the redshift range of interest (6 ^ z ^ 10), the rates of the two extreme models are similar, 
indicating that our conclusions are robust: GRB 080913 and GRB 090423 have Population II progenitors, and are both 
long bursts resulting from collapsars. 



